The Effect of Balancing Selection on Population Differentiation: A Study with HLA Genes.

Balancing selection is defined as a class of selective regimes that maintain polymorphism above what is expected under neutrality. Theory predicts that balancing selection reduces population differentiation, as measured by FST. However, balancing selection regimes in which different sets of alleles are maintained in different populations could increase population differentiation. To tackle the connection between balancing selection and population differentiation, we investigated population differentiation at the HLA genes, which constitute the most striking example of balancing selection in humans. We found that population differentiation of single nucleotide polymorphisms (SNPs) at the HLA genes is on average lower than that of SNPs in other genomic regions. We show that these results require using a computation that accounts for the dependence of FST on allele frequencies. However, in pairs of closely related populations, where genome-wide differentiation is low, differentiation at HLA is higher than in other genomic regions. Such increased population differentiation at HLA genes for recently diverged population pairs was reproduced in simulations of overdominant selection, as long as the fitness of the homozygotes differs between the diverging populations. The results give insight into a possible "divergent overdominance" mechanism for the nature of balancing selection on HLA genes across human populations

Evolutionary implications of a high selfing rate in the freshwater snail Lymnaea truncatula.

Self-compatible hermaphroditic organisms that mix self-fertilization and outcrossing are of great interest for investigating the evolution of mating systems. We investigate the evolution of selfing in Lymnaea truncatula, a self-compatible hermaphroditic freshwater snail. We first analyze the consequences of selfing in terms of genetic variability within and among populations and then investigate how these consequences along with the species ecology (harshness of the habitat and parasitism) might govern the evolution of selfing. Snails from 13 localities (classified as temporary or permanent depending on their water availability) were sampled in western Switzerland and genotyped for seven microsatellite loci. F(IS) (estimated on adults) and progeny array analyses (on hatchlings) provided similar selfing rate estimates of 80%. Populations presented a low polymorphism and were highly differentiated (F(ST) = 0.58). Although the reproductive assurance hypothesis would predict higher selfing rate in temporary populations, no difference in selfing level was observed between temporary and permanent populations. However, allelic richness and gene diversity declined in temporary habitats, presumably reflecting drift. Infection levels varied but were not simply related to either estimated population selfing rate or to differences in heterozygosity. These findings and the similar selfing rates estimated for hatchlings and adults suggest that within-population inbreeding depression is low in L. truncatula

Do riverine barriers, history or introgression shape the genetic structuring of a common shrew (Sorex araneus) population?

The common shrew (Sorex araneus) is subdivided into numerous chromosome races. The Valais and Cordon chromosome races meet and hybridize at a mountain river in Les Houches (French Alps). Significant genetic structuring was recently reported among populations found on the Valais side of this hybrid zone. In this paper, a phylogenetic analysis and partial Mantel tests are used to investigate the patterns and causes of this structuring. A total of 185 shrews were trapped at 12 localities. All individuals were typed for nine microsatellite loci. Although several mountain rivers are found in the study area, riverine barriers do not have a significant influence on gene flow. Partial Mantel tests show that our result is caused by the influence of the hybrid zone with the Cordon race. The geographical patterns of this structuring are discussed in the context of the contact zone, which appears to extend up to a group of two rivers. The glacier they originate from is known to have cut the Arve valley as recently as 1818. The recent history of this glacier, its moraine and possibly rivers, may therefore be linked to the history of this hybrid zone

Continental-scale patterns of pathogen prevalence: a case study on the corncrake

Pathogen infections can represent a substantial threat to wild populations, especially those already limited in size. To determine how much variation in the pathogens observed among fragmented populations is caused by ecological factors, one needs to examine systems where host genetic diversity is consistent among the populations, thus controlling for any potentially confounding genetic effects. Here, we report geographic variation in haemosporidian infection among European populations of corncrake. This species now occurs in fragmented populations, but there is little genetic structure and equally high levels of genetic diversity among these populations. We observed a longitudinal gradient of prevalence from western to Eastern Europe negatively correlated with national agricultural yield, but positively correlated with corncrake census population sizes when only the most widespread lineage is considered. This likely reveals a possible impact of local agriculture intensity, which reduced host population densities in Western Europe and, potentially, insect vector abundance, thus reducing the transmission of pathogens. We conclude that in the corncrake system, where metapopulation dynamics resulted in variations in local census population sizes, but not in the genetic impoverishment of these populations, anthropogenic activity has led to a reduction in host populations and pathogen prevalence

Genetic basis of adaptation in Arabidopsis thaliana: Local adaptation at the seed dormancy QTL DOG1

Local adaptation provides an opportunity to study the genetic basis of adaptation and investigate the allelic architecture of adaptive genes. We study DELAY OF GERMINATION 1 (DOG1), a gene controlling natural variation in seed dormancy in Arabidopsis thaliana and investigate evolution of dormancy in 41 populations distributed in four regions separated by natural barriers. Using F ST and Q ST comparisons, we compare variation at DOG1 with neutral markers and quantitative variation in seed dormancy. Patterns of genetic differentiation among populations suggest that the gene DOG1 contributes to local adaptation. Although Q ST for seed dormancy is not different from F ST for neutral markers, a correlation with variation in summer precipitation supports that seed dormancy is adaptive. We characterize dormancy variation in several F 2-populations and show that a series of functionally distinct alleles segregate at the DOG1 locus. Theoretical models have shown that the number and effect of alleles segregatin at quantitative trait loci (QTL) have important consequences for adaptation. Our results provide support to models postulating a large number of alleles at quantitative trait loci involved in adaptation. © 2012 The Author(s). Evolution © 2012 The Society for the Study of Evolution.Peer Reviewe

Mapping Bias Overestimates Reference Allele Frequencies at the HLA Genes in the 1000 Genomes Project Phase I Data.

Next-generation sequencing (NGS) technologies have become the standard for data generation in studies of population genomics, as the 1000 Genomes Project (1000G). However, these techniques are known to be problematic when applied to highly polymorphic genomic regions, such as the human leukocyte antigen (HLA) genes. Because accurate genotype calls and allele frequency estimations are crucial to population genomics analyses, it is important to assess the reliability of NGS data. Here, we evaluate the reliability of genotype calls and allele frequency estimates of the single-nucleotide polymorphisms (SNPs) reported by 1000G (phase I) at five HLA genes (HLA-A, -B, -C, -DRB1, and -DQB1). We take advantage of the availability of HLA Sanger sequencing of 930 of the 1092 1000G samples and use this as a gold standard to benchmark the 1000G data. We document that 18.6% of SNP genotype calls in HLA genes are incorrect and that allele frequencies are estimated with an error greater than ±0.1 at approximately 25% of the SNPs in HLA genes. We found a bias toward overestimation of reference allele frequency for the 1000G data, indicating mapping bias is an important cause of error in frequency estimation in this dataset. We provide a list of sites that have poor allele frequency estimates and discuss the outcomes of including those sites in different kinds of analyses. Because the HLA region is the most polymorphic in the human genome, our results provide insights into the challenges of using of NGS data at other genomic regions of high diversity

A Unified Characterization of Population Structure and Relatedness.

Many population genetic activities, ranging from evolutionary studies to association mapping, to forensic identification, rely on appropriate estimates of population structure or relatedness. All applications require recognition that quantities with an underlying meaning of allelic dependence are not defined in an absolute sense, but instead are made "relative to" some set of alleles other than the target set. The 1984 Weir and Cockerham [Formula: see text] estimate made explicit that the reference set of alleles was across populations, whereas standard kinship estimates do not make the reference explicit. Weir and Cockerham stated that their [Formula: see text] estimates were for independent populations, and standard kinship estimates have an implicit assumption that pairs of individuals in a study sample, other than the target pair, are unrelated or are not inbred. However, populations lose independence when there is migration between them, and dependencies between pairs of individuals in a population exist for more than one target pair. We have therefore recast our treatments of population structure, relatedness, and inbreeding to make explicit that the parameters of interest involve the differences in degrees of allelic dependence between the target and the reference sets of alleles, and so can be negative. We take the reference set to be the population from which study individuals have been sampled. We provide simple moment estimates of these parameters, phrased in terms of allelic matching within and between individuals for relatedness and inbreeding, or within and between populations for population structure. A multi-level hierarchy of alleles within individuals, alleles between individuals within populations, and alleles between populations, allows a unified treatment of relatedness and population structure. We expect our new measures to have a wide range of applications, but we note that their estimates are sensitive to rare or private variants: some population-characterization applications suggest exploiting those sensitivities, whereas estimation of relatedness may best use all genetic markers without filtering on minor allele frequency

Cirsium species show disparity in patterns of genetic variation at their range-edge, despite similar patterns of reproduction and isolation

Genetic variation was assessed across the UK geographical range of Cirsium acaule and Cirsium heterophyllum. A decline in genetic diversity and increase in population divergence approaching the range edge of these species was predicted based on parallel declines in population density and seed production reported seperately. Patterns were compared with UK populations of the widespread Cirsium arvense.Populations were sampled along a latitudinal transect in the UK and genetic variation assessed using microsatellite markers. Cirsium acaule shows strong isolation by distance, a significant decline in diversity and an increase in divergence among range-edge populations. Geographical structure is also evident in C. arvense, whereas no such patterns are seen in C.heterophyllum. There is a major disparity between patterns of genetic variation in C. acaule and C. heterophyllum despite very similar patterns in seed production and population isolation in these species. This suggests it may be misleading to make assumptions about the geographical structure of genetic variation within species based solely on the present-day reproduction and distribution of populations